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ABSTRACT 

We report the discovery of a luminous ultra-soft X-ray excess in a radio-loud narrow-line Seyfert 1 galaxy, 
RX J1633+4718, from archival ROSAT observations. The thermal temperature of this emission, when fitted 
with a blackbody, is as low as 32.5^q eV. This is in remarkable contrast to the canonical temperatures of 
~0. 1-0.2 keV found hitherto for the soft X-ray excess in active galactic nuclei (AGN), and is interestingly 
close to the maximum temperature predicted for a postulated accretion disc in this object. If this emission is 
indeed blackbody in nature, the derived luminosity (3.5^ J x 10 44 ergs" 1 ) infers a compact emitting area with 
a size (~5x 10 12 cm or 0.33 AU in radius) that is comparable to several times the Schwarzschild radius of a 
black hole at the mass estimated for this AGN (~3x 1O 6 M ). In fact, this ultra-steep X-ray emission can be 
well fitted as the (Compton scattered) Wien tail of the multi-temperature blackbody emission from an optically 
thick accretion disc, whose parameters inferred (black hole mass and accretion rate) are in good agreement with 
independent estimates using optical emission line spectrum. We thus consider this feature as a signature of the 
long-sought X-ray radiation directly from a disc around a super-massive black hole, presenting observational 
evidence for a black hole accretion disc in AGN. Future observations with better data quality, together with 
improved independent measurements of the black hole mass, may constrain the spin of the black hole. 
Subject headings: accretion disks - galaxies: active - galaxies: jets - galaxies: Seyfert - X-rays: galaxies - 
galaxies: individual (RX J1633+4718) 



1. INTRODUCTION 

Accretion of matter onto a black hole (BH) via an optically 
thick disc (Shakura & Sunyaev 1973; Novikov & Thorne 
1973) is widely accepted as the engine to power active galac- 
tic nuclei (AGN) and black hole X-ray binaries (BHXBs) 
in a bright state. In observation, a search of evidence for 
BH accretion discs has long been pursued, and one impor- 
tant signature to look for is the predicted characteristic multi- 
temperature blackbody radiation. For BHXBs (BH masses 
M~ 10M©) the bulk of the disc blackbody emission, which 
has a maximum temperature fcr max ~0.5-l keV, falls within 
the readily observable X-ray bandpass, and the disc model has 
proven to be prevailing (see e.g. Done et al. 2007a, for a re- 
cent review). In AGNs habouring super-massive black holes 
(SMBH) with typical M~ 1O 6_9 M and hence much lower 
?max, however, the situation is much less clear. The main dif- 
ficulty lies in that the bulk of the disc emission falls within 
the far-to-extreme UV regime, which is not observable due 
to Galactic absorption. While it is generally thought that the 
optical/UV "big blue bump" of AGN is likely associated with 
disc emission (e.g. Shields 1978; Malkan 1983), some key is- 
sues have not been settled yet, however (e.g. Laor et al. 1997). 

A long-standing puzzle is that, although the (Compton scat- 
tered) Wien tail of the disc blackbody emission is expected 
to plunge into the soft X-ray band (e.g. Ross et al. 1992; 
Shimura & Takahara 1995a), convincing evidence remains il- 
lusive. Indeed, a soft X-ray excess is commonly detected in 
AGNs; however, their thermal temperatures, which appear to 
fall within a narrow range {kT ~0. 1-0.2 keV), are too high 
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to conform the disc model prediction (e.g. Gierliriski & Done 
2004), even with the effects of Compton scattering taken into 
account (Ross et al. 1992). Neither the inferred luminosities 
agree with the model values self-consistently (e.g. Miniutti 
et al. 2009). Furthermore, the independence of temperature 
on M found recently argues against the disc blackbody ori- 
gin (e.g. Bianchi et al. 2009). In general, the soft X-ray ex- 
cess is suggested to have different origins from disc emission 
(e.g. Done et al. 2007b). Although there were limited pre- 
vious attempts to link the soft X-ray excess with disc emis- 
sion in a few AGNs, e.g. Mkn766 (Molendi et al. 1993) and 
RE J 1034+396 (Puchnarewtcz et al. 2001), whereby M and 
accretion rates were constrained from the X-ray (and opti- 
cal/UV) data, the evidence is far from convincing. 

In this paper, we report on an ultra-soft X-ray excess de- 
tected in an AGN RX J1633+4718, whose thermal tempera- 
ture {kT ~30 e V) is among the lowest ever detected. We show 
that this emission is well consistent with blackbody radiation 
from an optically thick accretion disc around a SMBH. We 
assume a cosmology with Hq= 70 km s" 1 Mpc" 1 , Q,m = 0.3, 
and f?A = 0.7. Errors are quoted at the 90% confidence level 
unless mentioned otherwise. 

2. RXJ1633+4718 

RXJ1633+4718 (SDSS J163323.58+471859.0) was first 
identified with a narrow-line Seyfert 1 (NLS 1 ) galaxy at a red- 
shift z=0.1 16 in optical identification of the ROSAT All-Sky 
Survey (RASS) sources (Moran et al. 1996; Wisotzki & Bade 
1997). The AGN resides in one of a pair of galaxies, with the 
other a starburst galaxy separated by 4" (Bade et al. 1995). 
It is associated with a variable, invert spectrum radio source 
(Neumann et al. 1994) with an unresolved VLBI core at milli- 
arcsec resolution (Doi et al. 2007). It is among a very radio- 
loud NLS1 AGN sample studied by Yuan et al. (2008). This 
rare type of AGNs exhibit properties characteristic of blazars, 
i.e. having relativistic jets directed close to the line-of-sight 
(Zhou et al. 2007; Yuan et al. 2008), which has been con- 
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firmed in several of the objects by recent 7-ray observations 
with Fermi (Abdo et al. 2009). 

Its optical spectrum, acquired in the Sloan Digital Sky Sur- 
vey (SDSS) and analysed in Yuan et al. (2008, see their 
Figure 1), has broad components of the H/3 and Ha lines 
with a width FWHM = 909 ±43 km s" 1 and a luminosity 
3.14(1 1.7) x 10 41 ergs" 1 , respectively. The expected "thermal" 
continuum luminosity at 5100 estimated from the H/3 lumi- 
nosity is AL^gioo = 2.63 x 10 43 ergs" 1 using the relations from 
Greene & Ho (2005). The "thermal" bolometric luminosity 
(i.e. free from the non-thermal jet emission) can then be es- 
timated as Lboi=£ALA5ioo =2.6(±0.5) x lO^ergs" 1 , adopting 
k = 10 ±2 (la) for radio-quiet AGNs (Richards et al. 2006). 
The BH mass, estimated from the broad Balmer linewidth and 
the luminosity of either the broad line or the "thermal" con- 
tinuum lies within M~ 2.0-3.5xl0 6 M Q using various com- 
monly adopted formalisms (e.g. Vestergaard & Peterson 2006; 
Collin et al. 2006; Greene & Ho 2007). We thus adopt a nom- 
inal value M= 3 x 10 6 M©, with a systematic uncertainty of 
~ 0.3 dex (Vestergaard & Peterson 2006). The estimated Ed- 
dington ratio is then L/LEdd=0.69^Q 35, using the above "ther- 
mal" bolometric luminosity. 

3. ULTRA-SOFT X-RAY EXCESS EMISSION 

J1633+4718 was observed with the ROSAT PSPC-b on axis 
with a 3732 s exposure on July 24th, 1993, and also detected 
during the RASS (PSPC-c) with a weighted exposure of 909 s. 
The reduction and a brief analysis of the ROSAT data was pre- 
sented in Yuan et al. (2008). To summarise, there are ~ 976 
net source counts in the pointed observation and 185 counts 
in the RASS. The 0. 1-2.4 keV spectrum can be well modeled 
with two components, an ultra-soft X-ray component domi- 
nating energies below 0.4 keV and a flat power-law (photon 
index T ~ 1 .37) dominating energies above 0.4 ke V; the latter 
was interpreted as inverse Compton emission from a relativis- 
tic jet (Yuan et al. 2008). In this paper we focus on the ultra- 
soft X-ray component and revisit the spectrum by performing 
a more detailed analysis. 

We use the pointed observation spectrum (Figure 1) to de- 
rive spectral parameters, given its much higher data signal-to- 
noise ratio (S/N) than that of the RASS data. XSPEC (v. 12.5) 
is used to perform spectral fit. Neutral absorption with a H I 
column density A^h is always added in spectral models con- 
cerned. The Galactic A^ al =1.79 x 10 20 cm" 2 on the line-of- 
sight, measured from the LAB Survey (Kalberla et al. 2005). 
The results are summarised in Table 1 . We first try a simple 
power-law model. The fit is unacceptable for fixing A^H=A^ al 
(X 2 /d.o.f.=67/30), underestimating severely the fluxes at both 
the low and high energies. Setting A^h free still yields a poor 
fit, since the same systematic structures remain in the residu- 
als and, even worse, resulting in little or no absorption. Fitting 
the higher-energy (0.4-2. 4keV) spectrum only (fixing Nu= 
A^ al ) yields a good fit and a flat T = 1.63(±0.4), while extrap- 
olating this power-law to lower energies reveals a prominent 
soft X-ray excess, as shown in Figure 1 (middle panel). 

Therefore an additional emission component is added to ac- 
count for the soft X-ray excess, that is modeled by commonly 
adopted models including blackbody, power-law, thin plasma 
and thermal bremstrahlung. In general, all these models yield 
nearly the same excellent fits and the systematic structures in 
the residuals vanish. However, they cannot be distinguished 
based on fitting statistics. For most of the models, and partic- 
ularly blackbody, the fitted A^h values are in good agreement 
with A^ al , albeit large uncertainties. We thus fix the absorp- 



tion Nh=Nu s1 when deriving spectral parameters (Table 1). As 
an example, the best-fit power-law plus blackbody model and 
the residuals are shown in Figure 1 . We note that the addition 
of a second, relatively steep power-law continuum as in nor- 
mal Seyferts with fixed V =2-3 is not required, indicating that 
it is negligible. 

Interestingly, the fitted (rest frame) thermal temperatures 
are surprisingly low, e.g. kT^ = 32.5+g[jeV for a (redshifted) 
blackbody. Such an exceptionally low temperature is in re- 
markable contrast to the "canonical" ^0. 1-0.2 keV found for 
AGNs (Gierlihski & Done 2004; Haba et al. 2008; Bianchi et 
al. 2009), and is perhaps the lowest with confident measure- 
ment known so far 4 (to our knowledge). In fact, the ultra- 
softness of this component is evident, as it dominates the 
emission below 0.4 keV only. The overall blackbody lumi- 
nosity is high, Lbb = 3.5^j'5 x 10 44 ergs _I , whose confidence 
contours vs. £7bb are shown in Figure 2. In fact, the black- 
body luminosity is comparable to the "thermal" bolometric 
luminosity [~ 2.8(±0.6) x 10 44 ergs _1 ], which is remarkable 
considering the large uncertainty in the bolometric correction. 
Hence we assume this ultra-soft X-ray emission as blackbody 
in nature. 

The spectrum measured in the RASS has an almost iden- 
tical spectral shape as that of the pointed observation. The 
best-fit parameters are in excellent agreement with those ob- 
tained from the pointed observation, with T = 1 .47+j] ^ and 
the blackbody temperature 30^jo eV (90% confidence level), 
though the errors are large given the small source counts. The 
overall blackbody luminosity is found to be 3.4 x 10 44 ergs _1 . 
There appears to be no detectable variation in either the spec- 
tral shape or the luminosity between the RASS and pointed 
observation. 

4. EMISSION FROM AN ACCRETION DISC? 
4.1. Size and temperature of X-ray emitting region 

Assuming blackbody radiation, the size of the emitting area 
can be inferred as A = L\,b/(aT^ h ) (a the Stefan-Boltzmann 
constant). For an assumed geometry of a sphere or a face-on 
disc, a radius ~ 5(or 7)x 10 12 cm (~0. 3-0.5 AU) is inferred, 
which is extremely compact for its high luminosity. In fact, 
these radii are merely several to ten times the Schwarzschild 
radius of a BH at the mass estimated for RX J 1633+47 18 (R s = 
2GM/c 2 sa 9 x 10 11 cm). An immediate conclusion is that the 
soft X-ray emitting region must be extremely compact, e.g. 
most likely in the close vicinity of the central BH. 

In the standard BH accretion disc model the local effective 
temperature at a radius R is 
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(Frank et al. 1992; Kato et al. 2008), where R m is the radius of 
the inner boundary (assumed to be the last marginally stable 
orbit) of the disc, and m = M /Mndd the scaled accretion rate, 
with the Eddington rate defined as Msdd = ^EddA^c 2 ), e.g. 
1.39 x 10 18 M/M© gs" 1 assuming rj = 0.1. For a non-rotating 

4 Lower kT in a few objects were reported by Urry et al. (1989) based on 
Einstein IPC+MPC data; however, they were not confirmed in latter observa- 
tions, e.g. in Mkn766 (Molendi et al. 1993). 
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BH = 3R S ), from outer part inwards T reaches its maxi- 
mum 

r ^ = 1L5 (io^) 1/4 " 1/4(eV) ' (2) 

at R = (49/36) x 3R S , and decreases again inwards. For 
M=3 x 10 6 M Q and assuming m as L/L E dd (=0.69), we find 
/tr max = 25.3 eV, that is close to the fitted r bb = 32.5^;g eV. To 
our knowledge, this is by far the most convincing case where 
the soft X-ray excess temperature is well compatible with the 
predicted maximum disc temperature. 

4.2. Accretion disc model fit 
4.2.1. Standard accretion disc model and parameters 

We fit the spectrum with the commonly adopted multi- 
colour disc (MCD, diskbb in XSPEC) model (Mitsuda et al. 
1984; Makishima et al. 1986). This model has been widely 
used in fitting the X-ray spectra of BHXBs, and found to 
provide reasonably good descriptions to the observed spec- 
tra at the high/soft state (e.g. Dotani et al. 1997; Kubota et 
al. 1998). With an approximation to the boundary condi- 
tion, the MCD model has a monotonic temperature profile, 
T e g(R) = Ti n (R/Ri n )~ 3 / 4 , where 7j n is the effective temperature 
at the disc inner radius This model should be sufficient 
for our purpose here, considering the moderate spectral S/N 
and resolution and the limited bandpass of the PSPC data. 

The emergent MCD spectrum is determined by two param- 
eters: the (apparent) radius of the inner boundary of the disc 
R[ n and the maximum (colour) temperature 7j' r which are re- 
lated to the disc luminosity via 

L = 4ttR'2<jT«. (3) 

and R' in vcos9 can be derived from spectral fitting, where 
9 is the disc inclination angle. Since the MCD model 
does not incorporate the redshift effect, which is small for 
RX J1633+4718 at z=0.1 16 though, T{ a and R' m in the object's 
rest frame are obtained by multiplying the fitted values with 
(1+z) and 1/(1 +z) (can be shown by using Eq. 3), respec- 
tively. The true inner boundary radius R m is related to the ap- 
parent radius R' m via R m = £R' in , where £ = ^/3jl{6 /7) 3 « 0.41 
is introduced to account for the inner boundary condition that 
is neglected in the MCD model (Kubota et al. 1998; Kato et 
al. 2008). 

The disc inclination 9 is expected to be a small angle given 
the blazar-like property of RX J1633+4718, i.e. generally < 
10° assuming disc-jet perpendicularity. In fact, for a wider 
range of 9, say 9 < 30° (as expected for type 1 AGNs), the 
dependence of \/co&9 on 9 is weak. We thus expect \/co&9 ps 
1, and hence /?i n vcos0ss7?i n . 

To relate the observed spectrum to the true disc parameters 
corrections should be made to account for the spectral hard- 
ening effect, which gives rise to a shifted Wien spectrum to a 
higher colour temperature due to incoherent Compton scatter- 
ing (Czerny & Elvis 1987; Wandel & Petrosian 1988; Ross et 
al. 1992; Shimura & Takahara 1995a). The true effective r max 
is then 

Tmax = 7] n , (^) 
K 

where k is the spectral hardening factor. Consequently, the 
true inner boundary radius 

Rm = k 2 £S£, (5) 



(Kato et al. 2008; Kubota et al. 1998). The value of k is some- 
what uncertain, and k ~ 1.7 ±0.2 was suggested (Shimura 
& Takahara 1995b) based mostly on consideration of BHXBs 
and 1 < k < 2.5 for SMBHs with wide ranges of M and m 
(Ross et al. 1992; Shimura & Takahara 1995a). Here we adopt 
1.7. 

4.2.2. Results of accretion disc model fit 

The spectrum is well fitted with the MCD plus a power- 
law model, with the fitted A^h consistent with A^j . Fixing 
N H =N^ yields the same excellent fit as using blackbody, 
giving £7^=39. 6!} 2 5 3 eV and R[ n =3 .0^ x 10 12 /\/cos6»cm in 
the object's rest frame (Table 1). Applying the corrections 
in Eqs. 4 and 5 with k ~ 1.7, we find a true maximum ef- 
fective temperature kT m . dx =23 .3 J ^\ eV and an inner disc ra- 
dius R in =3.6+{ 4 7 x 10 12 /\/cos7J cm. The confidence contours 
of kT mdx and R\ n are shown in Figure 3. 

Assuming a non-rotating BH and /?;„= 3R S , we can derive 
its mass from X-rays, M x , and in turn m using r max and Eq. 2. 
We find the best estimates M x = 4.1 x 10 6 /Vcos#M Q and 
m= 0.68/v^osTJ. The contours of confidence intervals for M 
and in, as derived from those of R m and T lwdx assuming 6> = 10°, 
are shown in Figure 3 (right panel). It is interesting to note that 
the X-ray spectrum provides relatively tight constraint on ac- 
cretion rate m. This is because the uncertainties of the derived 
^in (M x ) and T mdx , albeit large, are highly inversely coupled 
(left panel) and hence largely canceled mutually. 

These derived parameter values can be compared with those 
estimated from the optical spectrum, M and L/Ledd (assumed 
as in), and in turn R m and T max (using Eq. 2 and assum- 
ing R ln =3R s ). Their best estimates and uncertainty ranges 
are overplotted in Figure 3, as ±0.3 dex in M (Vestergaard 
& Peterson 2006) and ±2(1 cr) in the bolometric correction 
(Richards et al. 2006). Remarkably, the results derived from 
the X-ray and optical observations agree strikingly well, albeit 
the relatively large uncertainties in both the measurements. 
This conclusion is basically not affected by the exact value of 
k, provided 1 < k < 2.5. 

We note that the MCD model is only approximately correct 
since no general relativistic (GR) effects are taken into ac- 
count. We make further corrections to the derived disc param- 
eters due to the GR effects following the approach of Zhang 
et al. (1997). Eqs.4 and 5 are then T mdx = /"V" 1 ^) and 
R m = g ~°- 5 f 2 (K 2 ^R' ia ), where / and g are temperature and flux 
correction factors, respectively. For a non-rotating BH and 
9 w 0°, / w 0.85 and g w 0.80 (Zhang et al. 1997, their Ta- 
ble 1). Therefore the above kT nYdx and (M x ) values should 
be further multiplied by a factor of 1.17 and 0.81, respectively, 
i.e. A:r max =27eV,fl in ~2.9x 10 12 cm, M x ~ 3.3 x 10 6 M Q , and 
m^ 1.0. These parameter values are still marginally consis- 
tent with those derived from the optical spectrum (Figure 3). 

Therefore, we conclude that we might be seeing the direct 
X-ray emission from an accretion disc around the SMBH in 
RX J1633+4718, which has been long sought for AGNs. We 
note that, however, for a more realistic modeling, a relativis- 
tic disc model taking the GR effects fully into account should 
be used (such as the kerrbb model of Li et al. 2005). Be- 
sides, we did not consider the BH spin here, that is also plau- 
sible; e.g. for a spinning BH and a prograde disc, R m would 
extend further inwards within 3R S , i.e. ^i n <3iJ s , and r max in- 
creases, resulting in an increased M x and a reduced m. Fur- 
thermore, since the above found m>0.3, the effects of a slim 
disc (Abramowicz et al. 1988, Sadowski 2009, Li et al. 2010, 
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see Chapter 10 of Kato et al. 2008 for a review) may start to 
set in. We defer all these comprehensive treatments to future 
work, in a hope that it may shed light on the spin of the BH. 

5. DISCUSSION 
5.1. Compatibility with lower-frequency data 

Having established the accretion disc emission model from 
the soft X-ray spectrum, it would be interesting to check 
whether this model is compatible with the observed broad- 
band spectral energy distribution (SED) of RX J1633+4718. 
The unfolded PSPC spectrum and the fitted MCD plus power- 
law model (dotted line) are shown in Figure 4, along with the 
broad-band SED whose low-energy data are collected from 
the NASA/IPAC Extragalactic Database (NED). Interestingly, 
while the UV-to-infrared spectrum steepens towards lower 
frequencies, the Galex far-UV measurement (squares) agrees 
roughly with the MCD model prediction. The lack of an ob- 
vious big-blue-bump in optical-UV supports that it is likely 
shifted to higher frequencies. As for the optical-infrared 
bump, we interpret it as the high-energy tail of beamed syn- 
chrotron emission from the relativistic jets, in light of the 
blazar-like property of RX J1633+4718. We model this com- 
ponent with a parabolic function, a commonly used approx- 
imation to the synchrotron emission of blazars. The sum of 
this model and the above best-fit MCD model is fitted to the 
radio 5 GHz, near-IR and optical/UV data, excluding the far- 
IR measurements from IRAS and ISO (crosses; since they 
are most likely seriously contaminated, or even dominated 
by thermal emission of dust 5 from the companion starburst 
galaxy 4"away given the relatively poor spatial resolution of 
IRAS and ISO). As can be seen, such a synchrotron component 
(dashed line) can generally account for the SED at lower ener- 
gies. Therefore, the above accretion disc model, together with 
the other two emission components (synchrotron and inverse- 
Compton), is compatible with the observed broad-band SED 
of RX J1633+4718. Moreover, the lack of detectable varia- 
tion of the ultra-soft X-ray emission between the RASS and 
pointed observations also supports the accretion disc model, 
since the thermal emission from an accretion disc is generally 
stable over such a relatively short timescale. 

We note that the starburst companion must fall within the 
source extraction region on the PSPC detector plane due to 
its low spatial resolution. Starburst also produces soft X-rays, 
which are related to the far-IR luminosity as £p}^/£o.3-2keV = 
10 3-4 . For the companion starburst galaxy, an upper limit 
on Lpf R can be set by the IRAS measurements using /ioo M = 
1.12 Jy and/ 60A1 = 0.595 Jy (Sanders & Mirabel 1996), yield- 
ing Z|f R < L FIR = 3.0 x 10 n L Q (1.15 x lO^ergs" 1 ). This 
corresponds to an upper limit on the soft X-ray luminosity 
arising from the starburst Lf^ 3 _ 2keV < 1 x 10 41-42 ergs" 1 only. 
Such a contribution, if any, would be overwhelmed by the 
X-ray emission from the AGN given the observed luminos- 
ity L(j.3-2kev ~ 10 44 ergs _1 . Moreover, the thermal tempera- 
tures of X-ray gases of starburst fall within a strict range of 
0.3-0.7 keV, far higher than the fitted 50 eV for a thin plasma 
model (Table 1). We thus conclude that the ROSAT spectrum 
of RX J1633+4718 is essentially not affected by any potential 
contamination from the starburst companion. 

5.2. Why is RXJ1633+4718 different from other AGNs? 

5 The observed /(60/i)//(100/i)=0.53 is "infrared warm", typical of star- 
burst in the local universe (Soifer et al. 1987). 



An interesting question arises as to why RXJ1633+4718 
is so special and unique, in terms of previous unsuccessful 
attempts to ascribe the soft X-ray excess (with the canoni- 
cal temperatures of ~0. 1-0.2 keV) to the direct thermal emis- 
sion of the accretion discs in the vast majority of AGNs. 
This question will not be answered until the origin of the 
puzzling soft X-ray excess in many other AGNs is under- 
stood. Nonetheless, we discuss some possibilities here, albeit 
somewhat speculative. In any case, one important feature of 
RXJ1633+4718 is its dominance of the disc thermal emis- 
sion, which is a close analogue to the typical high/soft state 
in BHXBs (relatively weak corona; see above). This is differ- 
ent from many other AGNs in which the soft excess was in- 
vestigated, where the spectral state is characterised by strong 
corona emission. 

One of the viable explanations of the soft excess is Comp- 
tonisation of soft photons (e.g. from discs) off a region with 
electron temperature kT e ^0.1 -0.2keV and a large optical 
depth r ~ 20 (e.g. Czerny & Elvis 1987; Gierlihski & Done 
2004; Done et al. 2007b). In this scenario, RXJ1633+4718 
may differ from the other AGNs in its environment surround- 
ing the central engine, i.e. lacking this Comptonisation region 
and the disc being directly seen in the line-of-sight, which 
may not be the case for the other AGNs. For instance, the 
difference may be due to the presence of jets (and the face-on 
inclination of the disc) in RX J1633+4718. It may be pos- 
sible that in radio-quiet AGNs the Comptonisation region is 
formed by aborted jets (Henri & Petrucci 1997; Ghisellini et 
al. 2004) whose kinetic energy is converted to internal en- 
ergy to heat electrons to the required temperature. In the case 
of RX J1633+4718, the jets are successfully launched, how- 
ever, clearing the way to the central disc along the line-of- 
sight. It is interesting to note that in 3C273, a well-known 
AGN with relativistic jets, a relatively low blackbody temper- 
ature of the soft excess was also reported (~ 60 eV; Grandi & 
Palumbo 2004; although in that work the accretion disc origin 
was not tested in terms of matching the X-ray data with the 
predicted temperature and luminosity of the disc, as we did 
here). Whether or not this can be generalised to other radio- 
loud AGNs needs further investigations, though it should be 
reminded that the jets of radio-loud NLSls might be differ- 
ent from those of classical radio quasars, e.g. with reduced 
physical size and power (Yuan et al. 2008). 

An alternative model is disk reflection in which the ob- 
served soft X-ray excess can be explained as relativistically 
blurred line emission of the reflection component from a 
highly ionized inner disc, that somehow dominates the pri- 
mary emission (Ross & Fabian 1993; Ballantyne et al. 2001; 
Miniutti & Fabian 2004; Crummy et al. 2006). The postula- 
tion of this model is the presence of the primary X-ray emis- 
sion, presumably from a disc corona of high-temperature and 
low-optical depth. It may be possible that in RX J1633+4718 
the disc corona is significantly suppressed somehow, for in- 
stance, due to the formation of the jets (the energy goes into 
the jets instead of to the corona) or to the relatively high mass 
accretion rate. This would lead to a week reflection compo- 
nent in RX J 1633+47 18 and thus has no significant soft excess 
in the usual form, while the disc thermal emission stands out. 
In this line, as mentioned above, the Seyfert-like X-ray com- 
ponent in the ROSAT spectrum of RX J 1633+47 18 is too weak 
to be detected. 

Finally, for AGNs with high black hole masses (M > 
10 7 M Q ) and relatively low accretion rates, even if the 
other conditions (jets, corona, etc.) resemble those of 
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RX J1633+47 18, the thermal disc emission may have too low 
temperatures to be detectable above the low-energy cutoff of 
commonly used X-ray detectors (~ 0.1 keV). Future obser- 
vations of more AGNs similar to RX J1633+4718 may shed 
light on under what conditions the direct thermal emission of 
accretion discs can be detected in X-ray. 

6. CONCLUDING REMARKS 

We have shown that the ultra-soft X-ray excess discovered 
in RX J 1633+47 18 can be well and self-consistently described 
by thermal emission from an optically thick accretion disc 
around a SMBH. The derived parameters (M x , m, and T nmx ) 
agree with the independent estimates based on the optical 
spectrometric data. We thus consider this emission as a sig- 
nature of X-rays from an accretion disc around a SMBH, pre- 
senting possible evidence for BH accretion discs in AGNs. 
As such, RXJ1633+4718 is an AGN analog of BHXBs at 
the high/soft state in a sense that the accretion disc emission 
dominates the bolometric luminosity. More importantly, one 
may infer the spin of the BH, as succeeded in BHXBs (e.g. 
Zhang et al. 1997; McClintock et al. 2006). This may be 
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achieved with future improved observations, which, if suc- 
cessful, would provide a testbed for studying the link of BH 
spin and the formation of relativistic jets in AGNs. Finally, 
it is intriguing why RXJ1633+4718 is so unique among the 
many AGNs observed in the soft X-rays, which may be ad- 
dressed with future investigations. 
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TABLE 1 
Results of spectral fits 
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il Source model (as in XSPEC) modified by neutral absorption (wabs) with column 
density Nu: powl - power-law with photon index T; raymond - Reymond-Smith 
thin plasma emission (abundance fixed at solar); zbremss - redshifted thermal brem- 
strahlung; zbbdy - redshifted blackbody with effective temperature kT; diskbb - multi- 
colour disc model with the (rest frame) maximum temperature kT (see text) 

b r s : photon index of a second softer power-law 

c the Galactic /Vh value, fixed in fitting 



ill 
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Energy (keV) 

FIG. 1. — Upper panel: ROSATPSPC spectrum of RX J1633+4718 and the best-fit model of a power-law plus a blackbody (/V H fixed at Galactic /V* al ). Middle: 
data to model ratio, where a power-law model with fixed /V^ al is fitted to the 0.5-2.4 keV band and extrapolated down to 0. 1 keV. An overwhelming soft X-ray 
excess is evident. Lower: residuals of the best-fit model as in the upper panel. 
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FIG. 2. — Confidence contours (68%, 90%, 99% for two interesting parameters) for the fitted effective temperature and luminosity for the blackbody model 
component. 
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FIG. 3. — Confidence contours (solid; 68%, 90%, 99%) for the maximum effective temperature vs. the inner radius of the disc (left panel) and the BH mass 
versus scaled accretion rate (right panel; assuming a disc inclination 9 = 10°), derived from fitting the spectrum with the MCD plus a power-law model. The 
estimated values (asterisk) derived from the optical emission line spectrum and their uncertainty ranges are indicated (dashed; ±0.3dex in BH mass and ±2(lcr) 
in the bolometric correction factor are used). 
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FIG. 4. — Broad-band SED of RXJ1633+4718: the unfolded PSPC data (crosses) and the photometric measurements in the UV and optical (Galex and SDSS, 
squares) near-IR (2MASS, triangles), far-infrared (IRAS and ISO, crosses; arrows indicate upper limits) and radio (circles) bands. The PSF magnitudes are used 
whenever available. The models are the best-fit MCD+power-law model (dotted), synchrotron emission (dashed) and the sum (solid). 



